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ABSTRACT 

The nearby radio galaxy Centaurus A is poorly studied at high frequencies with conventional 
radio telescopes because of its very large angular size, but is one of a very few extragalactic 
objects to be detected and resolved by the Wilkinson Microwave Anisotropy Probe (WMAP). 
' We have used the five-year WMAP data for Cen A to constrain the high-frequency radio spec- 

C/3 , tra of the 10-degree giant lobes and to search for spectral changes as a function of position 

along the lobes. We show that the high-frequency radio spectra of the northern and southern 
giant lobes are significantly different: the spectrum of the southern lobe steepens monoton- 
, ically (and is steeper further from the active nucleus) whereas the spectrum of the northern 

J> ■ lobe remains consistent with a power law. The inferred differences in the northern and south- 

ed | ern giant lobes may be the result of real differences in their high-energy particle acceleration 

0^ ■ histories, perhaps due to the influence of the northern middle lobe, an intermediate-scale fea- 

' ture which has no detectable southern counterpart. In light of these results, we discuss the 

prospects for Fermi Gamma-ray Space Telescope detections of inverse-Compton emission 
OO ' from the giant lobes and the lobes' possible role in the production of the ultra-high energy 

cosmic rays (UHECR) detected by the Pierre Auger Observatory. We show that the possibil- 
ity of a Fermi detection depends sensitively on the physical conditions in the giant lobes, with 

■ the northern lobe more likely to be detected, and that any emission observed by Fermi is likely 

to be dominated by photons at the soft end of the Fermi energy band. On the other hand we 
argue that the estimated conditions in the giant lobes imply that UHECRs can be accelerated 

■ there, with a potentially detectable 7-ray signature at TeV energies. 
H ' 

Key words: radio continuum: galaxies - galaxies: jets - acceleration of particles - cosmic 
rays 



1 INTRODUCTION 

Centaurus A is the closest radio galaxy to us (we adopt D = 3.7 
Mpc, the average of 5 distance indicators in Ferrarese et al. 2007). 
Its proximity makes it one of the brightest extragalactic radio 
sources in the sky at low frequencies (only exceeded by Cygnus 
A: Baars et al. 1977), but also means that the outer 'giant' double 
lobes (throughout the paper we use the nomenclature adopted by 
Alvarez et al 2000) subtend an angle of ~ 10° on the sky, although 
their total physical size (~ 600 kpc in projection) is not unusually 
large for an radio galaxy of Cen A's luminosity. The large angular 
size of the lobes has prevented the type of spatially resolved, multi- 
frequency study of their spectral structure that is commonplace for 
more distant radio galaxies (e.g. Alexander & Leahy 1987). While 
detailed low-frequency maps of the giant lobes have been avail- 
able for many years (e.g. Cooper, Price & Cole 1965), published 
radio data from ground-based observations only exist up to 5 GHz 



(Junkes et al. 1993) and the spectral study of Alvarez et al. (2000), 
involving (at many frequencies) painstaking graphical integration 
of contour maps, was only able to determine overall spectra for the 
giant lobes, finding no evidence for deviation from a single power 
law between 408 MHz and 5 GHz. The low-frequency two-point 
spectral index maps of Combi & Romero (1997), however, do show 
some evidence for position-dependent spectral steepening, particu- 
larly towards the end of the southern giant lobe. 

Cen A is widely believed to be a restarting radio galaxy, in 
the sense that the inner lobes are the result of the current nuclear 
activity, while the giant outer lobes are the result of a previous 
outburst (Morganti et al. 1999). This picture is supported by the 
observation of hot thermal X-ray emission, apparently the result 
of strong shocks, surrounding both inner lobes (Kraft et al. 2003, 
2007; Croston et al., in prep.) which implies that they are propagat- 
ing supersonically into the intergalactic medium (IGM) of Cen A 
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and are disconnected from the giant lobes. However, the nature of 
the intermediate-scale northern middle lobe (NML: Morganti et al. 
1999) is not completely clear in this model. In standard spectral 
ageing models (e.g. Jaffe & Perola 1973), we might expect to see 
spectral steepening at high frequencies in the giant lobes; a mea- 
surement of spectral ageing gives a model-dependent constraint on 
the time since the last injection of high-energy electrons into these 
lobes. Such a constraint on spectral age could be compared with 
other estimates of the dynamical age of the radio source, and would 
therefore be of considerable interest, but the work of Alvarez et al. 
(2000) only sets upper limits on this quantity, since they do not see 
any deviation from a power-law spectrum. 

The giant lobes of Cen A are also interesting because they are 
predicted to be strong sources of inverse-Compton emission as the 
relativistic electrons in the lobes scatter cosmic microwave back- 
ground (CMB) photons to high energies; a detection of inverse- 
Compton emission from Cen A would constrain the magnetic field 
strength in the lobes of Fanaroff & Riley (1974) class I (here- 
after FRI) radio sources in general: we have little information on 
the magnetic field strengths in these low-power radio galaxies at 
present. However, X-ray emission from this process would be dis- 
tributed on similar scales to the giant lobes, making it hard to detect. 
Cooke, Lawrence & Perola (1978) claimed an early detection of the 
giant lobes using Ariel V, but Marshall & Clark (1981) argued that 
this was the result of point source contamination, placing a much 
lower upper limit on the flux from SAS 3 observations. At soft X- 
ray energies (e.g. Arp 1994) the situation is seriously confused by 
the presence of known thermal X-ray emission from the interstellar 
medium of the host galaxy, which more recently has been exten- 
sively studied with Chandra and XMM-Newton (e.g. Kraft et al. 
2003), and is also hard to study because the X-ray emission fills 
the field of view of modern soft X-ray imaging instruments such 
as ROSAT (Arp 1994), XMM, and ASCA (Isobe et al. 2001), pre- 
senting almost insuperable problems of background modelling and 
subtraction. However, the spectrum of the inverse-Compton emis- 
sion should be hard up to high energies (exactly how high de- 
pends on the model adopted for the electron energy spectrum, as 
we will discuss below). There are thus also interesting constraints 
from observations at MeV to GeV energies made with the Compton 
Gamma-Ray Observatory (e.g. Steinle et al. 1998; Sreekumar et al. 
1999), which do not detect the giant lobes but again set upper limits 
on their high-energy flux densities. More sensitive hard X-ray/7- 
ray observations exist with wide-field instruments like INTEGRAL 
(e.g. Rothschild et al. 2006) and the Swift Burst Alert Telescope 
(e.g. Markwardt et al. 2005) but as these are coded-aperture instru- 
ments they have limited sensitivity to extended emission (see e.g. 
Renaud et al. 2007). At present, therefore, there is no unambigu- 
ous detection of X-ray or 7-ray inverse-Compton emission from 
the giant lobes of Cen A. One of us has shown (Cheung 2007) that 
FermQ may have the sensitivity to detect inverse-Compton emis- 
sion off the CMB from the lobes of Cen A at energies from ~100 
MeV to 10 GeV. However, the details of this depend on modelling 
of the electron energy spectrum at high energies, which in turn de- 
pends on high-frequency radio data. 

Finally, Cen A's giant lobes are possible sources of ultra-high 
energy cosmic rays (UHECRs). Cen A's proximity means that all 
aspects of the active galaxy - central AGN, inner jets and lobes, 
and giant lobes - have long been considered as possible UHECR 
accelerators (see e.g. Cavallo 1978; Romero et al. 1996, and, more 



recently, Gureev & Troitsky 2008 and references therein). Interest 
in Cen A has been spurred by the remarkable discovery that 2 of 
the 27 UHECR events detected so far by the Pierre Auger Observa- 
tory (hereafter 'PAO' ; Abraham et al. 2007) appear to be arriving 
from the direction of the centre of Cen A, while at least 2 additional 
events may be associated with it (e.g. Gorbunov et al. 2008a; Wibig 
& Wolfendale 2007; Fargion 2008) due to the large angular extent 
of the giant radio lobes (Gorbunov et al. 2008b; Moskalenko et al. 
2008). Most scenarios discussed in the literature to date assume that 
UHECRs are produced near the supermassive black hole (SMBH) 
or in the inner jets (e.g. Cuoco & Hannestad 2008; Kachelriess, 
Ostapchenko & Tomas 2008), but an explanation in terms of the 
giant lobes has the advantage that it can easily account for the PAO 
events seen on larger scales. In order to investigate this quantita- 
tively we need information about the magnetic field strengths and 
the leptonic particle acceleration in these lobes, which can be pro- 
vided by a combination of high-frequency radio observations and 
inverse-Compton constraints or measurements. 

The Wilkinson Microwave Anisotropy Probe (WMAP) has ob- 
served the whole sky at frequencies around 20, 30, 40, 60 and 90 
GHz (known as K, Ka, Q, V and W bands respectively) with the aim 
of measuring structure in the CMB (e.g. Hinshaw et al. 2008). The 
currently available WMAP data represent 5 years of observations. 
Cen A is clearly detected, and spatially resolved, in the WMAP ob- 
servations at all frequencies (e.g. Page et al. 2007) and Israel et 
al. (2008) have recently presented WMAP-derived measurements 
of the flux density of the whole source, showing that there is clear 
steepening in the integrated spectrum at high frequencies. Thus the 
data are available to carry out a study of the variation of the radio 
spectrum as a function of position, to fit spectral ageing models to 
the large-scale lobes and investigate whether we can learn anything 
about the source dynamics, and to make predictions of the expected 
inverse-Compton emission from the giant lobes. 

In this paper we present the results of such a study. We first 
combine the 5-year WMAP data on Cen A with single-dish radio 
images at lower frequencies to make spatially resolved measure- 
ments of the radio spectra from 408 MHz to 90 GHz. We then dis- 
cuss the implications of the high-frequency detections for the dy- 
namics of Cen A, for possible inverse-Compton detections of the 
giant lobes at high energies, and for acceleration of UHECRs and 
their possible 7-ray emission signatures. 



2 DATA 

2.1 Radio data 

We obtained electronic versions of ground-based radio maps of 
Cen A to provide low-frequency counterparts to the five WMAP 
images. At 408 MHz we used the all-sky map of Haslam et al. 
(1982) which is available from the NCSA Astronomy Digital Im- 
age Librarjfl The 1.4-GHz and 4.75-GHz (hereafter, 5-GHz) data 
were taken with the Parkes telescope and kindly provided to us by 
Norbert Junkes; the 5-GHz data are those presented by Junkes et al. 
(1993). 

2.2 WMAP data and processing 

FITS maps of the 5-year WMAP data for Cen A at all 5 WMAP 
frequencies were kindly provided by N. Odegard. For each map, a 



Formerly known as the Gamma-ray Large Area Space Telescope, GLAST. 



2 See http://adil.ncsa.uiuc.edu/document/95.CH.01.01 . 
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Figure 1. Large-scale structure of Centaurus A with a resolution of ~ 0.83° . Contours are at 1, 2, 4 . . . times the base level specified for each map, and take 
no account of background, except for the 408-MHz map, from which a constant background of 43 Jy beam -1 has been subtracted. Top row, from left to 
right: 408 MHz (3 Jy beam" 1 ), 1.4 GHz (1 Jy beam" 1 ), 5 GHz (0.25 Jy beam" 1 ), 20 GHz (1.0 Jy beam -1 ). Bottom row, from left to right: 30 GHz (0.5 Jy 
beam -1 ), 40 GHz (1.0 Jy beam -1 ), 60 GHz (1.0 Jy beam -1 ), 90 GHz (2.0 Jy beam" 1 ). Circles in the bottom right-hand corner of each image indicate the 
beam size (diameter shows FWHM). The 20-GHz data have different noise characteristics from the other WMAP images because they have not been convolved 
with a Gaussian; instrumental noise on scales smaller than the effective beam is therefore visible. See the text (Section [2.2t for discussion of the convolution, 
effective resolution and beam area of these images. 



monopole defined outside of the WMAP Kp2 mask (see e.g. Hin- 
shaw et al. 2007) and the internal linear combination (ILC) CMB 
map were subtracted. The data were interpolated from HEALPix 
to a Galactic co-ordinate grid with 0.1° spacing. The images were 
then converted into a celestial co-ordinate frame using AIPS, and 
the units of the maps converted from mK to Jy beam -1 for ease of 
comparison with the ground-based radio maps, which have these 
units. For the unit conversion we used the beam areas tabulated by 
Page et al. (2003) and added a 'beam size' header to the FITS file 
that defined the beam to be a Gaussian of the same area. Although 
the WMAP beam is not a Gaussian (Page et al.), this approximation 
has no effect on our flux density measurements, since the two iden- 
tical values of beam area cancel; in addition, the Gaussian FWHM 
of the beam calculated in this way is similar to the true FWHM 
tabulated by Page et al. (2003), and so gives us a convenient way to 
characterize the resolution of the resulting maps. We were then able 
to use AIPS for simple analysis of the WMAP data and software 



based on the Funtools packag^f] for measurement of flux densities 
from matched regions of the ground-based and WMAP radio maps. 
The effective central frequencies of the WMAP bands for extended 
emission with a spectral index of 0.7 are 22.5, 32.7, 40.4, 60.1 and 
92.9 GHz (based on the bandpass characterization of Jarosik et al. 
2003) and we treat measured flux densities as being measured at 
these central frequencies in all subsequent analysis, although for 
simplicity we will refer to the different bands as the 20, 30, 40, 60 
and 90-GHz bands in what follows. 

In measuring flux densities from regions of Cen A with WMAP 
we need to take careful account of the foreground and background 
characteristics of the data. Firstly, the ILC subtraction described 
above only takes account of structure in the CMB on scales greater 
than 1°: this gives rise to real structure on small scales in the maps 
of this 10° source, which is a problem at higher frequencies and 
correspondingly high resolutions. Secondly, synchrotron emission 

3 Documented at http://hea-www.harvard.edu/saord/funtools/ ; see also 
Mandel, Murray & Roll (2001). 
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Figure 2. Map of two-point spectral index a between 1 .4 and 30 GHz in 
the giant lobes of Cen A. Only pixels where the total intensity exceeds 5 
times the off-source r.m.s. level are plotted. Contours are from the 1.4-GHz 
Parkes data at 1, 2, 4 ... Jy beam -1 . The resolution (0.83°) is indicated by 
the circle in the bottom left-hand corner. 

from the Milky Way provides a foreground which is both position 
and frequency-dependent, being strongest at low frequencies and 
towards the south of the Cen A field (e.g. de Oliviera-Costa et al. 
2008). 

For our analysis we convolved all the images with a Gaussian 
that gave us approximately the resolution of the 20-GHz WMAP 
image (0.83°) at all observing frequencies, with the exception of 
the 20-GHz data themselves and of the 408-MHz dataset, which 
already has a very similar intrinsic resolution. The convolution has 
two effects on our analysis: for the WMAP data, it helps to reduce 
any contamination from the unsubtracted real structure in the CMB; 
and, for all the datasets, it allows us to measure flux densities from 
extended regions with sufficient confidence that we are seeing the 
same structures at different frequencies. 

Fig.Q]shows a montage of the radio images used in this paper. 
Full-resolution WMAP images are presented by Israel et al. (2008). 



3 RESULTS 

3.1 Spectral index mapping 

To search for evidence for spectral differences in the large- 
scale structure of Cen A as a function of position we ini- 
tially constructed a map of the two-point spectral index (a = 
-log(S 1 /S 2 )/log(vi/v 2 )) between 1.4 and 30 GHz (Fig.©. 
Such a map can only be approximate given that the beams of the 
ground-based and WMAP data are only approximately matched, 
and it cannot take account of the varying position-dependent back- 
ground in the WMAP images, but it does provide an indication of 
whether the spectral steepening already observed by Israel et al. 
(2008) is position-dependent. The spectral index map gives a strong 
indication that the spectral index does indeed steepen as a function 
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Figure 3. Spectral extraction regions used for Table \T\ and Fig. [4] The 
greyscale shows the 20-GHz WMAP data. 

of distance from the 'nucleus' (which at this resolution includes 
the entire inner lobe and north middle lobe structures as well as the 
flat-spectrum core), confirming the lower-frequency (0.4—1.4 GHz) 
results of Combi & Romero (1997). 

3.2 Flux density measurements 

To investigate spectral changes in the giant lobes more quantita- 
tively we next divided the radio source into 5 regions encompassing 
all of the structure detected at 20 GHz. These regions are shown in 
Fig.[3]and for simplicity we number them 1-5 from north to south 
on the source. Regions 1 and 2 are respectively the outer and in- 
ner regions of the northern giant lobe, region 3 is dominated by 
the emission from the nucleus, inner lobes and north middle lobe, 
and regions 4 and 5 are respectively the inner and outer regions 
of the southern giant lobe. The regions are all substantially larger 
than the adopted WMAP beam sizes. For each region we defined 
two background regions lying adjacent to the source regions in the 
east-west direction. The background regions were used both to find 
a mean background flux density, which was subtracted from the to- 
tal in the source regions (thus removing contamination from our 
own Galaxy, whose strength varies with position on the maps) and 
to determine an off-source r.m.s. for each map, from which we es- 
timated the on-source noise. Flux measurements and background 
subtraction were carried out in the same way for all the ground- 
based and WMAP datasets. The measurements of flux density for 
each region and frequency, together with the error estimated from 
the background region r.m.s., are tabulated in Table [T] and plotted 
together in Fig. [4] 

We initially note that the total flux density listed in Table Q] 
the sum of all measured flux densities from the regions, agrees well 
with the values tabulated by Israel et al. (2008) (their table 1) within 
the uncertainties. Our flux densities are generally slightly lower, 
and the discrepancy is greatest at the higher frequencies where we 
have not attempted to measure fluxes for the whole source, but even 
here the differences do not exceed the errors quoted by Israel et al. 



The giant lobes of Centaurus A 5 



Table 1. Flux densities of the regions of Cen A as a function of frequency. The final column gives the sum of all measured quantities in the preceding columns. 
Blanks indicate combinations of frequency and region where significant emission was not visible by eye. 



Frequency Flux density (Jy) 



(GHz) 


Region 1 


Region 2 


Region 3 


Region 4 


Region 5 


Total 


0.408 


191.3 ± 10.8 


165.3 ± 12.0 


1074.0 ±18.8 


366.3 ±26.7 


238.5 ±30.7 


2035.2 ± 47.6 


1.4 


91.2 ±2.1 


74.6 ±3.0 


545.0 ± 1.3 


204.7 ± 4.3 


111.2 ± 6.0 


1026.7 ±8.3 


4.75 


30.3 ±0.3 


26.3 ±0.9 


273.5 ± 1.0 


104.3 ± 0.9 


47.0 ± 1.7 


481.4 ±2.4 


22.5 


5.7±0.5 


4.9 ±0.3 


71.3 ± 0.6 


19.6 ± 0.6 


6.3 ± 0.7 


107.9 ± 1.2 


32.7 


4.2 ±0.4 


3.4 ±0.3 


56.4 ±0.4 


13.5 ± 0.4 


3.9 ± 0.5 


81.3 ±0.9 


40.4 


3.5 ±0.5 


2.8 ±0.4 


48.8 ± 0.5 


10.2 ± 0.3 


3.0 ±0.5 


68.3 ± 1.0 


60.1 






37.0 ±0.8 


5.2 ± 0.6 


2.2 ± 0.8 


44.3 ± 1.3 


92.9 






28.1 ± 1.9 






28.1 ± 1.9 
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Figure 4. Flux density X frequency (uf v ) for the regions of Table [T] (see 
Fig. [3} as a function of frequency. 1 Jy Hz is 10 — 23 erg s _1 cm -2 . 



(larger than ours because they do not restrict themselves to regions 
where emission is clearly detected). The agreement between the 
two independent analyses strengthens the case that the flux density 
measurements are correct. We confirm the spectral steepening of 
the integrated source spectrum detected by Israel et al. 

Examination of the SEDs for individual regions (Fig.[4]l shows 
clear differences between the northern giant lobe (regions 1 and 2), 
the central region (3), and the southern giant lobe (4 and 5). The 



central region shows no evidence for spectral steepening except be- 
tween 5 and 20 GHz. At WMAP frequencies the spectral index a is 
very consistently about 0.65, close to the lower-frequency spectral 
index of 0.70 derived for the inner lobes by Alvarez et al. (2000). 
Given that the central region (region 3) in our analysis includes 
the flat-spectrum core (which is strongly variable, e.g. Israel et al. 
2008), the inner lobes, the north middle lobe and some extended 
lobe emission, it is not surprising that its spectrum is poorly fitted 
by a single power-law or curved spectrum. Probably it consists of 
(at least) a flat-spectrum component and a component with a spec- 
tral turnover between 5 and 20 GHz. 

The spectra of the southern giant lobe regions (regions 4 and 
5) remain relatively flat up to 5 GHz, but show clear steepening at 
high frequencies. Moreover, the spectrum of region 5 (the southern 
half of the southern giant lobe) is systematically steeper than that of 
region 4 (the northern half), even at low frequencies (as indicated 
by our spectral index mapping and also seen by Combi & Romero 
1997). This difference is statistically significant: considering just 
the spectral indices between 1.4 and 5 GHz, 04 = 0.55 ± 0.02 
while Q5 = 0.71 ±0.05. The systematic increase in spectral curva- 
ture as a function of frequency, seen particularly in region 4 (north- 
ern south giant lobe), is characteristic of spectral ageing models, 
and we explore fits of such models to the data below. 

On the other hand, the northern giant lobe regions (regions 1 
and 2) show no significant spectral differences. Their spectra both 
start to steepen between 1.4 and 5 GHz and then remain roughly 
constant (at around a — 1.0) at higher frequencies. These spectra 
are clearly different from those of the southern giant lobe. It does 
not seem possible to explain this difference in terms of inadequate 
background subtraction; the Galactic background that we subtract 
off is a power law, so that it is hard to see how any combination 
of partially subtracted background and source could produce the 
differences in northern and southern regions. The statistical differ- 
ences (given our adopted errors) between the northern and southern 
giant lobes are highly significant both at low and high frequencies. 
We explore the reasons for this difference below (Section[5]l. 

3.3 High-resolution imaging of the NML 

Finally we note that in the 90-GHz WMAP image, the only one to 
have the required resolution to distinguish between the inner lobes 
and the NML, there is a clear extension of the central point source 
(Fig.O which is exactly coincident with the NML. Its flux den- 
sity at 90 GHz is around 1.5 Jy which implies a comparatively flat 
spectrum (a ~ 0.6) between 5 and 90 GHz, measuring flux from 
matched regions of the full-resolution 5-GHz map of Junkes et al. 
(1993) and the 90-GHz WMAP image. This strongly suggests that 
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Figure 5. WMAP and radio images of the central parts of Cen A, show- 
ing the elongated structure that we associate with the northern middle 
lobe. Contours are from the full-resolution WMAP 90-GHz image, at 
0.75 X (1, 2, 4 . . .) Jy beam -1 ) overlaid on the 5-GHz radio map of Junkes 
et al. (1993). The area of the WMAP beam at this frequency corresponds to 
an effective resolution around 15 arcmin. The compact feature to the SW of 
the inner lobes/core of Cen A, undetected at 90 GHz, is a background radio 
galaxy. 



particle acceleration is either ongoing there or has recently ceased, 
although given that we only have two spectral data points we can- 
not rule out the possibility that the spectrum is flatter simply due to 
a higher magnetic field strength (cf. Katz-Stone, Rudnick & Ander- 
son 1993). We discuss the implications for the nature of the NML 
in Section|4] 



4 THE NATURE OF THE NML 

As discussed in Section Q] the northern middle lobe (NML: e.g. 
Morganti et al. 1999) is a feature at the base of the northern gi- 
ant lobe, peaking at around 0.5° (30 kpc in projection) from the 
active nucleus, that has no counterpart on the southern side. High- 
resolution radio data (e.g. Junkes et al. 1993; Fig. [5jl show a high- 
surface-brightness feature, centred at around RA = 13 h 27 m Dec 
= —42° 40', which extends smoothly northwards into the larger- 
scale lobe, with the magnetic field vectors oriented along the long 
axis of the structure (Junkes et al. 1993). At higher resolutions still 
the brightest part of the NML has been imaged by Morganti et al. 
(1999), who see a jet-like structure (the 'large-scale jet') connect- 
ing the NML and the northern inner lobe. Although high-resolution 
imaging of the northern inner lobe (e.g. Clarke et al. 1992) shows 
a sharply bounded radio structure with no evidence for continued 
outflow, it is possible (as argued by Morganti et al.) that the exis- 
tence of the large-scale jet implies that there is continued energy 
supply to the NML. This interpretation would be consistent with 
the flat spectrum between 5 and 90 GHz of the base of the NML 
region (see above). It would also be consistent with the detection of 
X-ray emission apparently embedded in the lobes (Feigelson et al. 
1981; Morganti et al. 1999) which is shown by XMM observations 



to be highly overpressured thermal material, plausibly requiring an 
interaction with an ongoing outflow (Kraft et al. 2008, submitted to 
ApJ.). The NML is a highly complex region associated with optical 
filaments and kinematically disturbed neutral hydrogen (Morganti 
et al. 1999; Oosterloo & Morganti 2005) and any interpretation is 
difficult, particularly as we are not aware of analogous structures 
in other restarting sources (Section[T}. However, our data support 
the idea that the NML is, or at least was until recently, fed by an 
active jet. Assuming an equipartition magnetic fielcQ, an 'injection 
index' (i.e. the synchrotron spectral index of the low-energy elec- 
tron population) of a = 0.5 and a Jaffe & Perola (1973: hereafter 
J-P) ageing spectrum, the data require particle acceleration to have 
ceased no more than 7 x 10 6 years ago. 



5 SPECTRA AND SPECTRAL AGEING IN THE GIANT 
LOBES 

As shown above, the spectra of the giant lobes of Cen A present a 
mixed picture. The southern giant lobe has a steep high-frequency 
spectrum and (within the limitations of our data) the spectra appear 
to steepen systematically both with frequency and with distance 
from the nucleus (Figs [2] @. In the northern giant lobe the spec- 
tra seem very uniform as a function of position and do not steepen 
significantly with frequency above 5 GHz. What causes these dif- 
ferences in the two giant lobes? 

We begin by noting that, given the large integration regions, 
the southern giant lobe spectra are remarkably well fitted by stan- 
dard J-P aged synchrotron spectra with an injection index of a — 
0.5 (determined from the low-frequency spectral indexfl We es- 
timate an equipartition field strength in the integration regions 
around 1.3 /iG (treating them as cylinders in the plane of the 
sky), which implies that inverse-Compton losses dominate, since 
the CMB has an energy density equivalent to a magnetic field with 
a strength of 3.3 /xG. If we assume that the field strength has the 
equipartition values and that the CMB is the dominant photon field, 
the spectral fits (see Fig. [6] for an example) imply spectral ages of 
(2.4 ± 0.1) x 10 7 years for region 4 and (2.9 ± 0.1) x 10 7 years 
for region 5 (errors are la statistical errors for one interesting pa- 
rameter only: systematic errors are dominated by the choice of in- 
jection index and ageing field strength). By contrast, regions 1 and 
2 (corresponding to the northern giant lobe) are rather poorly fit- 
ted by such models, since their low-frequency spectra are too steep 
and their high-frequency spectra too flat. They are better fitted with 
standard 'continuous injection' broken power-law models with an 
injection index of 0.5 steepening to 1.0 at high frequencies (Heav- 
ens & Meisenheimer 1986). If a J-P model is fitted the derived ages 
are, unsurprisingly, of the same order of magnitude but, implausi- 
bly, region 2 appears older than region 1 (ages (2.1 ±0.1) x 10 7 
and (2.9 ±0.1) x 10 7 years for regions 1 and 2 respectively, where 
we derive the age uncertainties by rescaling the errors on the data 
points so as to make the reduced \ 2 of the fit equal to 1). The ages 



4 Here and throughout the paper, except where otherwise stated, we assume 
an electron-positron plasma with no significant energetic contribution from 
protons. 

5 We note that the injection indices we assume are significantly flatter than 
the spectral index of the inner lobes measured by Alvarez et al.. This is not 
necessarily a problem, since 1 ) the equipartition magnetic field strength in 
the inner lobes is much higher and 2) there is no particular reason to believe 
that the particle acceleration process currently operating in the inner lobes 
is the same as that which operated in the giant lobes. 
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we derive are similar to the upper limits quoted by Alvarez et al. 
(2000), who used somewhat lower magnetic field strengths. 

One possible explanation for the spectral differences between 
the two giant lobes is that the northern lobe has undergone some 
particle injection event that the southern lobe has not. A natural 
explanation is that this is connected to the existence of the NML 
(Section[4]l: in this model the NML would be currently (or very re- 
cently) connected to the energy supply and injecting high-energy 
electrons into the base of the giant lobe. In this picture, the bright 
radio structure extending from the NML into the northern giant 
lobe on larger scales could be interpreted as continued outflow of 
recently accelerated electrons, although this would have to extend 
to 100-kpc scales to fully explain the peculiar spectra of lobe re- 
gions 1 and 2. The strong differences in the polarization structures 
of the giant lobes (Cooper et al. 1965; Junkes et al. 1993), and in 
particular the fact that the northern giant lobe is strongly polarized 
with the magnetic field direction aligned along the long axis of the 
NML, support this picture, as do differences in the fine structure of 
the radio emission in the northern and southern giant lobes (Feain 
et al. in prep). 

If this is the case, then the spectral ages derived for the south- 
ern giant lobe (where there is no evidence of a connection to the 
energy supply) allow us to suggest that the last acceleration of the 
electrons now dominating the synchrotron radiation in this region 
took place around 3 x 10 7 years ago (with the usual large system- 
atic uncertainties due to the assumption of equipartition). This does 
not seem unreasonable given that various estimates of the dynam- 
ical age of the inner lobes give a few x 10 6 years. The dynamical 
age of the giant lobes, if they have expanded at speeds of the order 
of the sound speed in the hot external medium, with kT ~ 0.35 
keV (Kraft et al. 2002), is expected to be almost an order of mag- 
nitude higher, but this is consistent if we assume that high-energy 
particles continued to be supplied to the southern giant lobe by a 
(now dissipated) jet throughout its lifetime. 



6 PROSPECTS FOR INVERSE-COMPTON 
DETECTIONS 

The mixed results on the spectra of the giant lobes imply ambigu- 
ous prospects for inverse-Compton detections at high energies. We 
estimate the inverse-Compton spectrum for an equipartition mag- 
netic field using the code of Hardcastle, Birkinshaw & Worrall 
(1998), taking the flux densities of the giant lobes to be the sum 
of the fluxes for regions 1 and 2 (north lobe) and 4 and 5 (south 
lobe) in TableQ] and treating the lobes as cylinders in the plane of 
the sky. We then use two reference models for the electron energy 
spectrum: a J-P spectrum with a spectral age of 3 x 10 7 years, as 
found for regions 4 and 5, for the south lobe, and a broken power 
law, which provides a better fit to regions 1 and 2, for the north 
lobe. 

The energetically dominant photon population, and therefore 
the population that dominates the inverse-Compton luminosity, is 
provided by the CMB. At the large distances of the Cen A giant 
lobes the energy density in photons provided by the host galaxy of 
Cen A is negligible, and we do not include it in our calculations; if 
any inverse-Compton emission arises from scattering of the galactic 
photons it will preferentially be seen at small distances from the nu- 
cleus and will be hard to resolve from the nuclear high-energy emis- 
sion. However, as Georganopoulos et al. (2008) have shown, the 
ubiquitous extragalactic background light (EBL) provides a second 
photon population that must be taken into account in high-energy 
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Figure 6. Flux densities for region 4 (Table[T} fitted with a standard Jaffe- 
Perola (1973) aged synchrotron spectrum assuming an injection index (low- 
frequency spectral index) of 0.5, an ageing 73-field of 1.3 fiG and a spectral 
age of 2.4 X 10 7 years. 

inverse-Compton calculations. The spectral energy distribution of 
the EBL is not well known (indeed, Georganopoulos et al. propose 
inverse-Compton observations as an additional constraint on its 
properties) and so we must adopt a model to make a prediction. We 
base our estimates of EBL-related inverse-Compton emission on 
the spectral energy distributiorjfl derived by Raue & Mazin (2008): 
this, while not a complete model, is intended to be consistent with 
all the existing direct and indirect limits. The model we use is 
significantly lower in photon energy density at all wavelengths 
than the upper range of the models adopted by Georganopoulos 
et al. (2008) and taken from Mazin & Raue (2007). Our inverse- 
Compton code includes an approximation to the Klein-Nishina cor- 
rection to the effective cross-section of the electron, since these ef- 
fects start to become significant when the highest-energy photons 
in the EBL model (which cuts off at A = 0.1 ^m) are scattered 
by the highest-energy electrons in the lobes, and also account for 
the kinematical upper limit on the energy of scattered photons (i.e. 
hv < ym e c 2 ). Since these calculations are only approximate, the 
results at the very highest energies should be taken as indicative 
only. 

The results of our calculations are shown in Fig. |7j For the 
south giant lobe, where the J-P spectrum is used, we see a sharp 
cutoff in the CMB inverse-Compton prediction between about 10 21 
and 10 22 Hz (4-40 MeV). This is as expected, since in the J-P spec- 
trum there are no electrons above ~ 3 x 10 11 eV (corresponding to 
7 = 6 x 10 5 ) and significant depletion of the spectrum below that, 
and the gain in energy for the inverse-Compton process is of order 
7 2 . Thus we do not predict a Fermi detection for equipartition field 
strengths for the south giant lobe, either via CMB photons or via 
the EBL (Fig.[7J upper right panel). 

On the other hand, for the broken power-law spectrum used to 

6 See also http://www.desy.de/~mraue/ebl/ . 
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Figure 7. Inverse-Compton predictions for the giant lobes of Cen A. Radio data points are the sum of our measurements from the ground-based radio maps 
and WMAP data for the regions corresponding to the northern (left) and southern (right) giant lobes as described in the text. High-energy data points are the 
same for all panels: the soft-X-ray limit of Marshall & Clark (1981) is converted to a flux density at 2 keV (~ 5 X 10 17 Hz) on the assumption that a = 0.5. 
The 0.05 - 30 MeV (~ 10 19 - 5 X 10 22 Hz) points are the reported Compton Gamma-Ray Observatory (CGRO) OSSE and COMPTEL detections of Cen A 
(Steinle et al. 1998), obtained from the NASA Extragalactic Database (NED), while the bow-tie plot centred on 100 MeV is the EGRET spectrum reported 
by Sreekumar et al. (1999), and the point at 190 GeV (5 X 10 25 Hz) represents the upper limit reported by Aharonian et al. (2005) from HESS observations, 
converted to a flux density by assuming power-law emission with photon index 3. We treat the high-energy data points as strict upper limits on the emission 
from the giant lobes. We also plot an estimate of the point-source Fermi sensitivity after one year (dashed line). The solid line shows the predicted synchrotron 
emission for the specified electron spectrum, the dotted line shows the inverse-Compton prediction for CMB photons, and the dot-dashed line shows the 
prediction for scattering of the EBL. For the northern giant lobe we plot the synchrotron and inverse-Compton spectra corresponding to a broken power-law 
electron energy spectrum with two possible high-energy cutoffs, 4 X 10 11 eV (upper panel) and 1 X 10 11 eV (lower panel), as discussed in the text. For the 
southern giant lobe (right-hand panel) we plot the synchrotron spectrum and inverse-Compton predictions for a J-P model as described in the text, with the 
upper panel showing the predictions for an equipartition magnetic field and the lower panel those for a field strength a factor 3 below equipartition. 
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model the north giant lobe, the maximum energy of the electrons is 
unconstrained. If we were arbitrarily to assume a sufficiently large 
value ( > 10 13 eV) we would obtain a spectrum from scattering of 
the CMB that does not cut off even in the Fermi band (10 22 - 10 25 
Hz). However, the predicted spectrum would be inconsistent with 
the 30 MeV - 2 GeV EGRET detection of Cen A (Sreekumar et al. 
1999; Hartman et al. 1999). The derived EGRET position is coinci- 
dent with the nucleus so we can treat this measured spectrum, par- 
ticularly at low energies, as a strict upper limit on the spectrum of 
either of the giant lobes, since the angular resolution of EGRET is 
~ 1° or greater at these energies (see e.g. fig. 1 of Funk et al. 2008). 
If we adopt a lower value of i? m ax, ~ 4 x 10 eV, we obtain an 
inverse-Compton spectrum that does not violate the EGRET limit 
but still allows a detection of scattered CMB photonfl of emission 
from this process at the extreme soft end of the Fermi band. This 
i? max is the one used in Fig. [7] (upper left panel). Although still 
lower values of i5 m ax are possible, it cannot be reduced too much 
before it starts to predict an (unobserved) high-frequency cutoff in 
the synchrotron spectrum. However, values of ~ 10 11 eV are per- 
mitted by the data, and this predicts no emission in the Fermi band 
(Fig. [7] lower left panel). Again, neither of these models predicts 
strong emission from inverse-Compton scattering of the EBL in 
the Fermi band if it is close to the level of our adopted model, al- 
though it would be detected in the most sensitive part of the band if 
its level were close to the upper value used by Georganopoulos et 
al. (2008). 

The inverse-Compton calculation in Fig.UJassumes equiparti- 
tion, while in general we find in the lobes of FRIIs that the observed 
inverse-Compton emission exceeds the equipartition prediction by 
a factor of a few (e.g. Croston et al. 2005). Reducing the magnetic 
field strength both increases the number density of electrons (since 
the synchrotron emissivity is fixed) and increases their maximum 
energy (since the observed cutoff frequency is fixed). However, the 
existing observations, particularly the soft X-ray limit of Marshall 
& Clark (1981), constrain the magnetic field to be no more than a 
factor ~ 3 below equipartition. If we reduce the magnetic field by 
this factor, then the inverse-Compton prediction implies detectabil- 
ity at very soft energies of the scattered CMB even for the south- 
ern giant lobe (Fig. UJ lower right panel), although we emphasize 
that Fermi-LKT observations are challenging in this energy range 
due to the decreased effective area and an increase in the back- 
ground diffuse emission. In addition we expect to see quite signifi- 
cant, spectrally distinct emission from inverse-Compton scattering 
of the EBL. This prediction may already marginally conflict with 
the CGRO data, and more extreme electron-density excesses over 
the equipartition values are certainly ruled out both by the Mar- 
shall & Clark limit and the CGRO spectrum, but we can conclude 
that even for the southern giant lobe there is a region of parameter 
space with B < B cq where a Fermi detection is possible. 

The results to be expected from Fermi therefore depend sen- 
sitively on the magnetic field strength in the lobes and on the un- 
known high-energy electron spectral cutoff in the north giant lobe. 
Fig. UJ illustrates some of the possibilities. If the magnetic field 



' We use Fermi 1-year 5<j point-source sensitivities taken from 
http://www-glast.slac.stanford.edu/software/IS/glast_lat_performance.htm: 
these are appropriate fiducial values in the absence of observational 
constraints on Fermi's real performance for an extended source. Below 
100 MeV, the background is expected to be high and our extrapolation is 
based on the sensitivity curve for the "inner galaxy" presented in Funk et 
al. (2008). 



strength is significantly below equipartition, or if the high-energy 
cutoff is high enough in the north lobe, then we expect a Fermi 
detection of scattered CMB photons for our adopted sensitivities, 
although in all cases the Fermi emission must be very soft to avoid 
violating the limits imposed by EGRET between 30 MeV and 1 
GeV. If Fermi does detect the giant lobes, then we will measure 
the magnetic field strength, and the sensitivities are such that we 
should be able to observe in the Fermi spectrum the difference 
in the high-energy electron spectra of the two giant lobes implied 
by the WMAP data: the resulting constraints on the magnetic field 
strength will then allow a measurement of, or limit on, the spectral 
energy density of the background of photons at higher energies (e.g. 
from the EBL) in the lobes. If extended inverse-Compton emission 
is not detected, then Fermi will set field-dependent limits on the 
high-energy cutoff in electron energies in the northern giant lobe 
and will improve significantly on the existing limits on magnetic 
field strength, but the problem of detecting inverse-Compton emis- 
sion from Cen A, and thus measuring the magnetic field strength in 
the giant lobes, will remain unsolved. 



7 UHECRS FROM THE GIANT LOBES AND THEIR 
POSSIBLE RADIATIVE SIGNATURES 

7.1 Introduction 

As mentioned in SectionQ] at least 4 of the 27 PAO UHECR events 
in the energy range E p ~ (0.6 — 0.9) x 10 20 eV (Abraham et al. 
2007, 2008) may be associated with Cen A (e.g. Moskalenko et 
al. 2008). Although evidence has been presented that the UHECRs 
include some heavier nuclei as well as protons (e.g. Unger et al. 
2007), the level of uncertainty in the composition does not affect 
our order of magnitude estimates if, as we assume hereafter, all the 
events are identified with protons. In this section we argue that the 
giant lobes of Cen A may confine such extremely relativistic pro- 
tons, and may indeed be the sites for their acceleration. Our spectral 
analysis of the giant lobes suggests that there is ongoing continuous 
particle (electron) acceleration in the northern part of the source, 
requiring, in turn, the presence of magnetic turbulence which can 
energize different cosmic ray species via stochastic particle-wave 
interactions. In addition, the parameters of the giant radio lobes of 
Cen A estimated in the previous sections allow us to consider the 
problem of acceleration of UHECRs in Cen A in a more quantita- 
tive way than has previously been possible. 

7.2 Acceleration of UHECRs 

The Larmor radius of a cosmic ray with energy E p = E20 10 20 
eV gyrating in a magnetic field with intensity B = B-e 10 -6 G is 
r-L = E p /eB ~ 100 E 2 qBzI kpc. Thus, the giant lobes of Cen A, 
with an estimated equipartition magnetic field B_s ~ 1 and a spa- 
tial scale R ~ 100 kpc, satisfy the standard Hillas (1984) criterion 
for a possible source of UHECRs, since < R for E20 < 1. 
The lobes of radio galaxies are highly magnetized, rarified cavi- 
ties inflated by jets in the surrounding medium, and are expected to 
contain a predominantly tangled magnetic field component with the 
maximum wavelength of the turbulent modes possibly as large as 
the scale of the system. Since the giant lobes are plausibly close to 
equipartition between magnetic field and particle energy density, as 
indicated by our analysis in Sections[5]and|6] the velocities of such 
modes are expected to be very close to the speed of light. These 
are excellent conditions for the efficient stochastic acceleration of 
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protons injected by the jets into, and confined within, the lobes, 
by means of resonant Fermi-type processes. The shortest possible 
timescale for such acceleration corresponds to Bohm-type diffu- 
sion (where the particle mean free path is comparable to the parti- 
cle gyro-radius: we assume this timescale in what follows), and is 
roughl)[f] t acc ~ 10 tl/ c. This implies that 

tacc ~ 3.5 E 20 BZ\ Myr (1) 

which is comfortably smaller than the age of the lobes estimated 
in the previous section, tiobc ~ 30 Myr. On the other hand, the 
characteristic timescale for the diffusive escape of particles from 
the lobes is of the same order or longer only if E20 < 1, namely 

tdifi ~ 3 R 2 /r h c ~ 0.9 R 2 100 B-e Myr (2) 

where flioo = R/WO kpc. The maximum energy of cosmic rays 
available in the acceleration scenario we are considering, given by 
the condition t acc ~ t^ff (since all the other loss timescales are 
much longer than t^g, see below), is E20 ~ 0.5 B-e R100 (thus 
satisfying the Hillas criterion, since the condition t acc ~ idiff re- 
duces exactly to tl ~ R for the assumed acceleration and diffu- 
sion timescales). For the estimated lobe parameters B-e ~ 1 and 
R100 ~ 1, this maximum energy is in very good agreement with 
the observed energies of the PAO events associated with Cen A (see 
above). 

We next consider whether Cen A is powerful enough to ac- 
count for the detected flux of UHECRs. We begin by estimating 
the energetics of the giant lobes assuming rough equipartition be- 
tween the magnetic field, relativistic electrons and (not necessar- 
ily relativistic) protons. For a cylindrical volume matching the gi- 
ant radio structure with V = tv R 2 h ~ 5.5 x 10 71 R-ioo ^600 
cm 3 (the radius R = R100 100 kpc and height h = heoo 600 
kpc), Etot ~ 3 17 B V ~ 6.6 x 10 58 Bl 6 R 2 00 heoo erg, where 
Ub = B 2 /8n ~4x 10~ 14 B- 6 erg cm -3 is the magnetic energy 
density. This leads to a very rough limit to the power of the jets 
supplying energy to the lobes 

Lj ~ Etot/2t\ ohe ~ 3.5 x 10 43 B 2 _ & R 2 00 heoo ergs -1 (3) 

This estimate should be considered as a lower limit, as it neglects 
any mechanical work done by the jets on the surrounding medium 
and the actual jet lifetime may also be substantially shorter than 
the age of the giant lobes; however, it is of the same order of mag- 
nitude as kinetic power estimates for the jets powering the inner 
lobes of Cen A (see Kraft et al. 2003) and with jet power estimates 
for more distant FRI radio galaxies (e.g. Laing et al. 2002). The 
UHECR flux detected by PAO from a point source with N associ- 
ated events above Eth = 60 EeV energies, assuming a power-law 
form of the cosmic ray spectrum with a energy index s — 2, is 
-Fuhe ~ 0.33 x 10 17 NE^o km" 2 yr" 1 eV~~\ when corrected 
for the relative exposure appropriate for Cen A's declination (e.g. 
Cuoco & Hannestad 2008). The particle energy index s is poorly 
known in the energy range we are considering; for simplicity, we 
assume hereafter the standard value of s = 2, consistent with the 

8 In general we can write r. acc = r^/r/c, where the factor lj accounts 
for less than ideal efficiency of particle acceleration; see, e.g., Aharonian 
et al. (2002). In the specific case of stochastic acceleration by magnetic 
turbulence (e.g. Stawarz & Petrosian 2008) the acceleration timescale (in 
the Bohm limit and with a turbulence spectrum W(k) going as k~ ) is 
given by t acc = (i) A /c)" 2 (C/o/f/ tur b)(r L /c) where v A is the Alfven 
speed and Uo/U tU ib is the ratio of energy densities in the unperturbed 
magnetic field and the turbulent magnetic component. Thus for relativistic 
strong turbulence (da c > ^turb ^ Uo), we expect the factor I/77 ~ 10. 



electron injection spectrum within the giant lobes of Cen A used 
in our analysis, which results in equal power stored per decade of 
cosmic ray energy. For N — 4, the monochromatic particle energy 
flux is then [E 2 F] V he ~ 0.66 x 10~ 12 erg cm -2 s" 1 and the 
monochromatic particle luminosity is thus 

Luhe ~ 4tt d 2 [£^]uhe ~ 10 39 erg s" 1 (4) 

Even if we extrapolate the cosmic ray spectrum down to the lowest 
energies with s — 2, the total cosmic ray luminosity is found to 
be Ltot ~ IOLuhe ~ 10 40 erg s _1 , and thus it constitutes a 
negligible fraction ( < 0.1%) of the power supplied by the jets to 
the giant lobes. 

7.3 Radiative signatures 

As discussed in the previous subsection, the timescales and ener- 
getics are both consistent with the hypothesis that UHECRs are 
accelerated continuously within the giant lobes of Cen A. In this 
scenario, the ultrarelativistic protons within the extended and mag- 
netized lobes will radiate, with the main loss processes being due 
to synchrotron radiation and proton-proton (p — p) interaction^- 
The former process is characterized by a timescale 

i syn ~ 1.4 x 10 6 E^o 1 B~_l Myr (5) 

resulting in a synchrotron continuum peaking at photon energies 
£ syn ~ 25 E 2 q B-a keV (Aharonian 2002). In the latter case, the 
characteristic timescale is 

i pp ~ 1.7 x 10 6 nZ\ Myr (6) 

where n t h = n-4 10 -4 cm -3 is the number density of the cold 
gas within the giant lobes, and 7-rays with energies e 7 < 10 EeV 
for the maximum proton energy E p ~ 10 20 eV will be generated 
(see e.g. Aharonian 2002). 

It follows immediately from the above results (eqs [2] [4] O 
that the synchrotron emission of the cosmic-ray protons produced 
within Cen A's giant lobes is negligible, since it is expected to peak 
around photon energies ~ 10 keV with a luminosity of the order of 

L syn ~ Z/uhe idiff Asyn ~ 10 33 R-ioo B- 6 ergs" 1 (7) 

This is orders of magnitude below the expected X-ray luminosity 
of the giant lobes' electrons due to inverse-Comptonization of the 
CMB photon field for B- 6 ~ 1 (~ 2 x 10 41 erg s" 1 ), and large 
increases in the magnetic field strength over the equipartition value 
would be required for the proton synchrotron emission to become 
dominant. 

9 UHECR will also interact with the photon field within the lobes, via 
both photo-meson production and Bethe-Heitler pair production. The key 
quantity for the former process is the number density of photons with en- 
ergies above the photo-meson production energy threshold, 300 MeV in 
the proton rest frame, which correspond to an observed photon energy 
e* = 3 E^Q meV, or a frequency u* = 7 X 10 11 Hz (Aharonian 
2002). However, for E20 ~ 1, the numerically dominant photon field in 
the Cen A lobes around and above u* (by a very large factor) is the CMB, 
and it is well known (e.g. Greisen 1966) that the loss timescale is > 10' 
years even for E20 ~ 1, and much larger for lower energies. Although 
the number density of synchrotron photons available for pair production 
for a given UHECR energy is higher (because of the lower energy thresh- 
old for pair production) CMB photons are still numerically dominant and 
we know that pair production is negligible compared to photo-meson pro- 
duction when CMB losses are considered. We therefore neglect these loss 
processes in what follows. 
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The expected 7-ray emission from p — p interactions might, 
however, be somewhat more promising in terms of allowing an in- 
dependent detection of the high-energy protons. We note first that 
the timescale for the diffusive escape of cosmic ray protons from 
the system becomes longer than the timescale for p — p interac- 
tions below energies E p ~ 100 R 2 00 B-6 n-4 TeV. Thus, all of 
the power channeled to such particles in the acceleration process 
during the lifetime of the source is released as 7-ray emission, be- 
low photon energies e 7 ~ 10 TeV, with an expected photon index 
T 7 = s (Kelner, Aharonian & Bugayov 2006). Unfortunately, the 
number density of thermal protons within the lobes of radio galax- 
ies is not known. The lack of observed internal depolarization in 
low-frequency radio observations (Cioffi & Jones 1980) has been 
used to set temperature-independent limits of the order of 10~ 4 
cm" 3 for extended components of other radio sources (e.g. Eilek 
et al. 1984; Spangler & Sakurai 1985), but this method depends 
on assumptions about the field geometry (Laing 1984) and in any 
case the required multi-frequency radio polarization observations 
of Cen A are not available to us. Isobe et al. (2001) claim an ASCA 
detection of soft X-ray emission coincident with the outer part of 
the giant lobes, with kT = 0.62 keV: assuming an abundance of 
0.1 solar, we find that their quoted soft X-ray (0.5-2.0 keV) flux 
corresponds to a density nth = 1-6 x 10" 4 cm -3 if it comes from 
a uniform thermal plasma internal to the lobe. This is a strict up- 
per limit on the gas density since emission may also be coming 
from gas external to the giant lobes (increasing the abundance re- 
duces the derived density). Similarly, the Marshall & Clark (1981) 
2-10 keV flux limit acts as a strict upper limit on the density of hot 
(kT — 2.5 keV) gas in the giant lobes, giving nth ~ 1.3 x 10" 4 
cm -3 if all the emission comes from a uniform thermal plasma with 
the same abundance. These limits on densities look comparable to 
the expected parameters of the intergalactic medium at hundred- 
kpc distances from the centre of the galaxy group hosting Cen A 
(i.e., from the Cen A host galaxy), and hence we assume n-4 ~ 1 
for illustration, although we note that gas with this density must be 
relatively cold if it is not to dominate the energy density of the giant 
lobes. 

With this value of n_4, and with the model cosmic ray spec- 
trum s — 2 anticipated above and normalized to the PAO flux, 
the expected monochromatic 7-ray luminosity from p-p interac- 
tions emitted over the whole volume of the extended giant lobes 
at photon energies < 10 TeV is 

L 1 ~ Luhe tiobc App ~ 10 34 n-4 erg s" 1 , (8) 

corresponding to a monochromatic (e.g., 10 TeV) flux energy den- 
sity of only [eS(e)]-, ~ L-jAnd 2 ~ 10" 17 erg cm" 2 s" 1 . This 
can be compared with the upper limit set on the TeV emission from 
the nucleus of Cen A by the existing 4.2-h HESS observation, F(> 
0.19TeV) < 5.68 x 10~ 12 phcm" 2 s" 1 (Aharonian et al. 2005; 
cf. Fig. [7}, which, with a 7-ray photon index T 7 = 2, corresponds 
to a monochromatic luminosity of [eL(e)] 0.19 TeV < 2.8 x 10 39 
erg s" 1 , or a flux of 1.7 x 10" 12 erg cm -2 s — 1 . Thus it seems that 
the expected 7-ray luminosity is much below the sensitivity limit of 
any current TeV instrument. Note that in addition the angular reso- 
lution of HESS in the 0.1 — 100 TeV photon energy range is < 0.1°, 
which is much better than the ~ 10° — 1° point spread function of 
EGRET (Section[6j in the photon energy range ~ 0.03 — 1 GeV 
Cen A therefore does not appear point-like to HESS, whereas it is 
effectively point-like at least at low energies with EGRET, so that 
the sensitivity limit of HESS appropriate for an extended source 
will be substantially higher than the quoted point source value. 



On the other hand, the estimated TeV flux of ~ 10" erg 
cm" 2 s" 1 should be considered as a very conservative lower limit. 
This is because we have assumed a rather flat cosmic ray spec- 
trum within the whole proton energy range. For any more realistic 
(steeper) spectrum for the accelerated cosmic rays at the highest 
energy range, the expected 7-ray flux increases. For example, as- 
suming s — 2.7 (as used by other authors) above some break en- 
ergy Ehr < Eth, and s — 2 below Ehr, we obtain the PAO flux at 
E p > -Eth with N — 4 events equal to 

[£ 2 F]uhe ~ 1.1 x 10~ 12 (£ p /60EeV)" ' 7 erg cm" 2 s" 1 (9) 

(see Cuoco & Hannestad 2008), and thus the monochromatic cos- 
mic ray luminosity Lett ~ 1.8 x 10 39 (E b[ /60 EeV)" 7 erg 
s" 1 at proton energies E p < Ehr- For example, with Ebr ~ 
10 ls eV (see in this context Kachelriess et al. 2008), one obtains 
Lcr ~ 3.2 x 10 40 erg s"\ which is more than one order of mag- 
nitude higher than our previous estimate, resulting in an expected 
monochromatic 7-ray energy flux at photon energies < 10 TeV of 
the order of [e5'(e)] 7 ~ 3 x 10~ 16 erg cm" 2 s" 1 . Yet another rea- 
son why this is a conservative lower limit is that the re-processing 
of higher-energy 7-rays, produced by the Ep > 100 TeV protons 
has not been taken into account in our analysis, and this may in- 
crease the expected 7-ray flux by a factor of at least few. Still, 
in order to reach the sensitivity limit of modern Cherenkov tele- 
scopes, the total power channelled to the cosmic rays with energies 
Ep < 100 TeV would have to be of the order of 10% of the total 
jet power, Lcr ~ 0.1 Lj, corresponding to a total energy stored in 
cosmic rays from the lobes of Ecu ~ 10 58 erg. 

The possibility of detecting the high-energy 7-ray emission 
from p — p interactions in the giant lobes of Cen A, even though 
not particularly supported by our analysis presented above, is an ex- 
citing one and worth some attention, since the large extent of these 
lobes means that they can be resolved at the highest TeV 7-ray 
energies using HESS. In this way the contribution from the inner 
regions of Cen A (including the inner, sub-pc to kpc scale jet, the 
host galaxy, and active nucleus) to the total 7-ray signal could be 
spatially resolved from the emission from the lobes. In this context 
it is important to note that our analysis of the leptonic emission gen- 
erated within the extended lobes (Section [6} indicates that we ex- 
pect no significant TeV radiation produced by ultrarelativistic elec- 
trons via inverse-Comptonization of either the CMB or EBL pho- 
ton fields, since the limits from synchrotron and inverse-Compton 
continuum indicate a cut-off in the lobes' electron spectrum around 
energies E e < 1 TeV, providing an effective upper limit on the en- 
ergy of any inverse-Compton scattered photons (Fig. [7J. Thus, if 
any TeV emission is detected from the giant lobes, it can be safely 
identified with being hadronic in origin. This is in contrast to the 
situation at the lower photon energy range probed by Fermi, which 
can also resolve the giant lobes from the central regions of Cen A, 
since the expected leptonic (inverse-Compton) emission produced 
from the lobes is expected to dominate at GeV photon energies over 
any hadronic emission. 

Obviously, UHECRs can in principle be accelerated in other 
components of the very complex Cen A system. Indeed, models for 
UHECR production on the very smallest scales (around the SMBH 
or in the inner jet) have been considered extensively in the liter- 
ature. These models are all based around the general idea that the 
accreting SMBHs in the centers of at least some active galaxies pro- 
vide a sufficient drop in potential for such an extreme acceleration. 
Indeed, the electromagnetic force associated with the black hole 
of mass Mbh rotating with angular momentum J in an external 
magnetic field (supported by the accreting matter) with intensity 
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B nuc , is AV ~ J B nuc /MBH c (Phinney 1983). The precise value 
for Bnuc is not known, but it is generally expected that the nuclear 
magnetic field energy density cannot be higher than (and may be 
equal to) the energy density of the accreting matter (e.g. Ghosh & 
Abramowicz 1997). Thus, for an object accreting at some fraction 
r\ of the Eddington rate 

B nuc ~ 6 x 10 4 r) 1/2 M~ 1/2 G (10) 

where M 8 = M BH /10 8 M©. The spin of SMBHs in general is 
also a very poorly known parameter. However, assuming the spin 
paradigm for efficient jet production (Blandford 1990), one can ex- 
pect SMBHs in elliptical-hosted radio galaxies to be spinning at the 
maximal rate, J ~ Jmax = Mbh cr g , where ra = G Mbh/c 2 is 
the Schwarzchild radius of the hole (see Sikora, Stawarz & Lasota 

2007) . These give AV ~ 10 18 rj 1/2 M 1/2 statvolt, and therefore 
the maximum available energy for a test particle with charge e ac- 
celerated in this potential drop is as high as 

£ max ~ 3 x 10 20 r, 1/2 M 8 1/2 eV (11) 

For Ms ~ 1 in Cen A (Marconi et al. 2006; Haring-Neumayer et 
al. 2006) and assuming an accretion efficiency of order r\ ~ 0.01 
(Evans et al. 2004), this leads to £ max ~3x 10 19 eV, in agree- 
ment with the energies of the associated PAO events. However, 
even in these models it is not clear where and how the involved 
electric circuit closes. It may close very close to the black hole, or 
along the edges of extended radio lobes (see e.g. Blandford 2008). 
In the former case UHECR protons may indeed be generated in 
the inner regions of Cen A (see, in this context, Kachelriess et al. 

2008) ; however, in the latter case, stochastic interactions with the 
magnetic turbulence taking place in the giant lobes, as discussed in 
this paper (and by others, e.g. Fraschetti & Melia 2008), or some 
other (though more speculative) mechanism involving reconnec- 
tion of the lobes' magnetic field (see in this context Benford & 
Protheroe 2008), may mediate the particle acceleration. Irrespec- 
tive of the detailed mechanism, only in the case of UHECR pro- 
duction in the giant lobes will there be a clear radiative signature 
potentially detectable at GeV-TeV photon energies, since any anal- 
ogous emission produced in the inner regions of the radio galaxy, 
even though possibly stronger due to the richer particle and photon 
environment (and therefore more significant p — p or proton-photon 
interactions), will always be confused with other 7-ray 'nuclear' 
components (e.g., Chiaberge et al. 2001), including direct emission 
from the inner jet, or the emission from a giant pair halo created 
by the re-processing of such emission (as discussed in the context 
of Cen A by Stawarz et al. 2006). The detection or non-detection 
of extended GeV-TeV emission from the giant lobes of Cen A can 
therefore in principle help to distinguish between these models. 



8 SUMMARY AND CONCLUSIONS 

We have presented the first high-frequency, spatially resolved study 
of the giant lobes of Centaurus A. We confirm recent findings (Is- 
rael et al. 2008) that the overall lobe spectra steepen at frequen- 
cies above 5 GHz, but in addition we have shown that the northern 
and southern giant lobes are significantly different at these high 
frequencies: the spectrum of the southern lobe steepens monoton- 
ically (and is steeper further from the active nucleus) whereas the 
spectrum of the northern lobe, after an initial steepening, remains 
consistent with a power law to within the limits of our data. We 
suggest that this indicates a real difference in the particle acceler- 



ation history of the northern and southern giant lobes, perhaps due 
to the influence of the poorly understood NML region. 

Few FRI sources have so far been studied at radio frequen- 
cies above a few tens of GHz, with more attention having been 
given to the hotspots and jets of powerful FRIIs (e.g. Hardcastle & 
Looney 2008 and references therein). The WMAP data for Cen A 
show that it is both possible and useful to study the dynamics of 
FRIs with high-frequency data. The Atacama Large Millimetre Ar- 
ray (ALMA), in combination with current and next-generation tele- 
scopes operating at lower frequencies, will allow spectral studies of 
large samples of well-studied FRI radio galaxies to be carried out 
relatively easily. It will be interesting to compare the results of such 
work with the results presented here for Cen A, and in particular to 
search for other radio sources that show strong spectral differences 
between the lobes at high frequencies. 

Our results also have implications for the high-energy astro- 
physics of Cen A. If the simplest interpretation of our spectral 
modelling is correct, then the southern giant lobe is a true relic 
with no ongoing particle acceleration, while the northern giant lobe 
may have an electron spectrum that extends to substantially higher 
energies. While Fermi detections of both lobes remain possible, 
the northern lobe is likely to be more clearly detectable, and all 
our Fermi predictions imply a relatively soft spectrum for inverse- 
Compton scattering of the CMB. Upper limits in the soft X-ray 
and from the EGRET detector on board the CGRO already rule out 
magnetic field strengths more than a factor of a few below equipar- 
tition: Fermi will either detect the giant lobes or substantially im- 
prove this limit. 

Finally, we have shown above that the giant lobes of Cen A are 
potential sites for acceleration of protons to the highest observed 
energies, E p ~ 10 20 eV, although we are required to adopt the 
most optimistic acceleration timescale, corresponding to the Bohm 
limit for the particle-MHD wave interactions, in order to be effi- 
cient enough to account for the PAO UHECR flux. We have also 
found that in this scenario the 7-ray radiative signatures accompa- 
nying the acceleration process, resulting from the interactions of 
ultrarelativistic protons with the thermal gas within the lobes, can 
possibly be observed by Cherenkov telescopes and/or by Fermi, 
thus permitting a test of this scenario. 
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